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The decay widths of the W and Z bosons that medi- 
ate the weak interaction are precisely predicted within 
the standard model (SM). At Born level the W width 
Tw and mass M\y are related through the precisely de- 
termined Fermi coupling constant, Gf- Beyond lead- 
ing order, higher-order electroweak (EW) and quantum- 
chromodynamic (QCD) corrections, Sew ~ —0.4% and 
Sqcd ~ 2.5% respectively, modify the relation such that 

IV = (1 + Sew + Sqcd) 0? ]■ The uncertainty 

on the SM prediction IV = 2091 ± 2 MeV is dominated 
by the uncertainty on M\y with smaller contributions 
from the uncertainties on the higher order corrections [2j • 
Uncertainties on SM parameters, such as the Higgs boson 
mass, affect this prediction very weakly, and so a mea- 
surement allows an unambiguous test of the SM that can 
also be used to constrain other SM parameters such as the 
V cs CKM matrix element 0] . The average of previously 
published direct measurements of Tw from pp collisions 
at the Tevatron [4( and e + e~ collisions at LEP-II 5] has 
a combined uncertainty of 2.7% with the most precise 
determination from a single experiment (ALEPH) hav- 
ing an uncertainty of 5.1%. The most precise indirect 
determination [6j of Tw from a measurement of the ratio 
^ = J[^Tz^*e+e-) nas an uncer tainty of 2%. 

This Letter presents the world's most precise direct de- 
termination of Tw from a single experiment. The anal- 
ysis uses W — » ev and W — > [iv data with respective 
integrated luminosities of 370 pb _1 and 330 pb _1 col- 
lected by the CDF II detector at the Fermilab Tevatron. 

Neutrinos are undetectable by the CDF II detector 
and hence the invariant mass of the W boson cannot 
be reconstructed. Tyy is therefore determined from a 
fit to the distribution of the W transverse mass Mt = 

J 2 (p T p T ~ Pt ' Vt)i wn ere pj, and p^ are the mea- 
sured transverse momentum of the charged lepton and 
the transverse momentum of the neutrino as inferred 
from the observed missing transverse energy, respectively. 
Events with Mt > Mw arise predominantly from a com- 
bination of non-zero W width and finite detector resolu- 
tion. The width component of the high- Mr line-shape 
falls off more slowly than the resolution component, al- 
lowing a precise Tw measurement from the Mt > Mw 
events even in the presence of systematic uncertainties 
on the resolution. 

The components of the CDF II detector relevant to this 
analysis are described briefly here: a more complete de- 
scription can be found elsewhere [7J . A silicon microstrip 
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detector [8] is used to measure the distance of closest ap- 
proach in the transverse plane, do, of charged particles to 
the beamline. The momenta of charged particles are mea- 
sured using a 96-layer open-cell drift chamber (COT) § 
inside a 1.4 T solenoid. Optimal pt resolution is ob- 
tained by constraining trajectories to originate from the 
beamline. Electromagnetic and hadronic calorimeters, 
arranged in a projective tower geometry, cover the pseu- 
dorapidity range < 3.64 [lcj . In the region |7y| < 1.0, a 
lead/scintillator electromagnetic calorimeter (CEM) [11( 
measures electron energies and proportional chambers 
embedded at the shower maximum provide further in- 
formation on shower shapes and positions. A system of 
drift chambers outside the calorimeters is used to identify 
muons in the region \rj\ < 1.0 [12j. 

W — ► ev candidate events are selected by a large trans- 
verse energy electron trigger, and the electron shower is 
required to have transverse energy E T > 25 GeV [13] in 
the CEM. The ratio of the energy measured in the CEM 
and the charged-track momentum measured in the COT, 
E/p, must satisfy 0.8 < E/p < 1.3. The ratio of energy 
deposited in the hadronic (HAD) and CEM calorimeter 
towers is required to satisfy ^hadZ-Bcem < 0.07. The 
electron shower must be contained within a fiducial re- 
gion of the CEM, away from calorimeter cell boundaries, 
and must have a typical electron lateral shower profile on 
the projection of the COT track in the CEM. Contamina- 
tion by Z — > e + e~ events is reduced by rejecting events 
with an additional high pt track of opposite sign charge 
pointing to an un-instrumented region of the calorimeter. 

W — > [iv candidate events are selected by a large pt 
muon trigger and are required to contain a COT track, 
well matched to a track segment in the muon chambers, 
with transverse momentum p T > 25 GeV. The energy 
deposited in the electromagnetic and hadronic calorime- 
ters must be consistent with the passage of a minimum- 
ionizing particle. Requirements on the track do and track 
fit x 2 arc imposed to reject background. Events consis- 
tent with cosmic rays or those with an additional high-py 
track consistent with Z — > [i + pT decays are removed. 

The existence of a neutrino is inferred from a trans- 
verse momentum imbalance. The missing transverse mo- 
mentum, p^ = — (pj, + it) , must satisfy p T > 25 GeV. 
The components of the recoil transverse energy vector 
u are defined as JV Ei sin #i(cos (pi, sin<^j), for calorime- 
ter towers i with |ry| < 3.64, excluding those traversed 
by and surrounding the charged lepton. u receives con- 
tributions from initial-state QCD radiation, underlying- 
event energy, final-state photon radiation, and overlap- 
ping pp interactions. To reduce backgrounds and improve 
transverse mass resolution, the recoil energy must satisfy 
u < 20 GeV. The W — > ev (W — > \iv) sample con- 
sists of 127432 (108808) candidate events in the range 
50 < M T < 200 GeV and 3436 (2619) in the high M T 
range of 90 < M T < 200 GeV. 

Since the W and Z bosons share a common production 
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mechanism and the momenta of Z bosons can be directly 
reconstructed from their decay products, Z — > £ + £~ de- 
cays are used to model the detector's response to W — > Iv 
events. Samples of Z — > e + e~ and Z — > candidates 
are selected by requiring two charged leptons, with the 
same requirements as the W lepton candidates, with the 
exception that the muon chamber track match require- 
ment is removed for one of the muons in the Z — > H + p~ 
pair. The di-lepton invariant mass is required to satisfy 
80 < M u < 100 GeV. Samples of 2909 Z -> e+e" and 
6271 Z — > events with recoil energy u < 20 GeV 

are used to determine the scale and resolution of the lep- 
ton energy and momentum measurements. A second set 
of Z — > e + e~ and Z — > /i + fx~ control samples is defined 
with the u cut replaced by a di-lepton transverse mo- 
mentum cut, < 50 GeV, in order to constrain the W 
boson's transverse momentum spectrum and to provide 
an empirical model of the recoil. 

The W boson Mt spectrum is modeled using a 
Monte-Carlo simulation. The CTEQ6M [Tj| parton 
distribution functions (PDFs) are used, and W bo- 
son invariant masses v§ are generated according to 
an energy-dependent Breit-Wigner distribution: a(s) ~ 

Higher order QCD ef- 
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M^/s) 2 + sT w /M^ 

fects are included by generating the W bosons with a 
Pt distribution from a NLO combined with rcsummation 
QCD calculation [l4| with the non-perturbative prescrip- 
tion of [l5j . Photon radiation from the charged lepton is 
simulated using a 0(a) matrix-element calculation [l6| . 
Corrections for EW box diagrams are applied from the 
calculation of [l7j |. 

The charged leptons and radiated photons are passed 
through a custom detector simulation that models in de- 
tail the energy loss due to ionization and bremsstrahlung. 
The simulation also includes a parametric model of the 
u measurement as a function of the boson p?, tuned on 
data as described below. The same kinematic and geo- 
metric cuts used to select candidate events in the data 
are applied to the simulation. The simulation produces 
Mt spectra for Tw values ranging from 1.0 to 3.0 GeV 
with M w fixed at 80.403 GeV 0. 

This measurement relies on the accurate modeling of 
the Mt distribution over a wide range. The most impor- 
tant sources of systematic uncertainty affecting the Mj 
shape arise from the charged-lepton energy and momen- 
tum scales and resolutions, the recoil modeling, and the 
presence of backgrounds. All systematic uncertainties are 
evaluated by varying parameters in the simulation and 
then fitting the resulting Mt spectra with the nominal 
spectra. Uncertainties have been calculated separately 
for the fit region M T ut < M T < 200 GeV for M T ut values 
ranging from 80 to 110 GeV. While the statistical uncer- 
tainty decreases as M T ut is lowered, the systematic un- 
certainty increases. A value of M™' = 90 GeV gives the 
smallest total uncertainty. Backgrounds are added to the 



simulation Mt spectra which are then normalized to the 
number of data events in the region 50 < Mt < 90 GeV. 

The COT momentum scale is determined from a fit 
to the Z — > fi + pT invariant mass distribution with the 
Z mass constrained to the world average value Q. A 
consistent COT momentum scale is obtained from fits 
to the invariant mass distributions of J/ip — > 
and T — > fi + ' pT events Q. The difference between 
the three determinations has a negligible effect on this 
analysis. The contribution to the uncertainty on Tw 
in the W — > [iv channel Ar^ arising from the 0.04% 
uncertainty in the COT momentum scale, due to the 
Z -> p^pT statistics, is AI^ = 17 MeV. 

By scaling the resolutions predicted by a GEANT [l8| 
simulation of the COT to match the observed di-muon 
invariant mass distribution in Z — ► p + p~ decays, we 
obtain a momentum resolution of (j(\/pt) = (5.4±0.2) x 
10~ 4 GeV -1 . A consistent a(l/pr) is also determined 
using the E/p distribution of the W — ► ev data. The 
combined uncertainties from the Z — > p + p~ and E/p 
fits for the COT resolution give AT 1 ^ = 26 MeV. 

The CEM energy scale and resolution are determined 
from fits to the Z — * e + e~ invariant-mass distribu- 
tion with the Z mass constrained to the world aver- 
age value [2j and to the E/p distribution of electrons 
in W ^ ev events, using the calibration of p described 
above. The scales determined from the two methods are 
consistent and are combined to form a weighted average 
with an uncertainty of 0.04%. The contribution to the 
uncertainty on Tw in the W — > ev channel Ar^ arising 
from this uncertainty is Ar^ = 17 MeV. 

The CEM resolution fits constrain the constant 
term k in the CEM resolution function a(E)/E = 
13.5%/V£r(GeV) © k The E/p and Z -► e+e" 

fit results differ by 1.6 standard deviations. They are 
combined and an uncertainty is assigned that spans both 
values, as well as the values obtained when the E/p fit 
region is varied, to give k = 1.1 ±0.4%. This uncertainty 
on the CEM resolution gives Ar^ = 31 MeV. 

Energy loss by electrons and photons in the solenoid 
coil and associated material prior to the CEM, as well as 
energy leakage into the hadronic calorimeter, are param- 
eterized based on the results of a GEANT simulation. In 
addition to these simulated sources of CEM non-linearity, 
an additional per-particle intrinsic non-linearity is deter- 
mined from the W — > ev and Z — > e + e~ data by fit- 
ting the E/p distribution in bins of Et- Its uncertainty 
gives Ar^ = 12 MeV, resulting in a total uncertainty 
of 21 MeV on Tyy from the uncertainties on the electron 
energy scale determination. Finally, uncertainties in the 
modeling of very low-energy photons and the amount of 
passive material prior to the COT give Ar^ = 13 MeV. 

The recoil transverse energy vector u is used to deter- 
mine Pt and hence Mt- Since u comes predominantly 
from initial-state QCD radiation, which is balanced by 
the W or Z boson pt, we form an empirical model by 
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parameterizing its response and resolution as a function 
of pip. The parameters of the model are varied according 
to the covariance matrices obtained in the fits to Z data. 
The resulting uncertainties on u from the recoil model 
give Ar^ (Ar^) = 54 (49) MeV. The uncertainty in 
the modeling of the distribution is determined by fit- 
ting the distribution in Z — > £ + £~ decays and results 
in a 7 MeV common uncertainty on Tw 

Several background processes can mimic the W signal. 
The process W — > tv — > £vvv has a signature similar 
to W — » Iv decays, but with lower M T - Z — > 
events, where only one lepton is identified, can be recon- 
structed as W candidates. These two backgrounds can 
be accurately determined from Monte-Carlo simulation. 
QCD multi-jet backgrounds arise when one jet mimics a 
charged lepton and another is mismeasured to produce 
an energy imbalance. Since the region with apparently 
low p T is enriched in QCD background, the background 
normalization is estimated from a fit to the p T distri- 
bution in events where the p T and low Mt cut are not 
applied. The background p T spectrum is taken from data 
events in which some of the charged-lepton identification 
cuts have been reversed, and the signal spectrum is taken 
from the simulation. A decay-in-flight (DIF) background 
to the W — » fiv signal arises when kaons or pions decay 
to /iv inside the COT, resulting in mismeasured muon 
momentum and a large \ 2 value between the COT hits 
assigned to the track and the fitted track trajectory. This 
background is estimated from a fit to the x 2 distribution. 
The background spectrum is taken from events with a 
large do and the signal from Z — > events, which 

have negligible background. The background fractions 
over the entire region 50 < Mr < 200 GeV are indicated 
in Figure 1. In the 90 < M T < 200 GeV fit region the 
total background fraction is 4.0 ± 0.2% (10.8 ± 0.3%) in 
the electron (muon) channel. Varying the background 
predictions within these overall normalization uncertain- 
ties, as well as varying their Mt shapes, causes varia- 
tions in the electron (muon) Tw of 32 (33) MeV. The 
backgrounds that contribute the most to the width un- 
certainty are the QCD multijet (DIF) backgrounds in the 
electron (muon) channel. 

We also investigate small systematic uncertainties due 
to PDFs, Mw, EW corrections, lepton identification, and 
acceptance. The uncertainty on Tw arising from PDFs is 
determined using the variations defined by the CTEQ6M 
PDF eigenvector basis [HI. The PDF error sets span a 
90% confidence interval so the resulting Tw shifts are 
divided by 1.6 to obtain la uncertainties [l!|, giving an 
uncertainty of 16 MeV in both channels. A systematic 
uncertainty of 12 MeV is added in quadrature to this to 
account for the effect of higher order QCD effects not 
implemented in the Monte-Carlo simulation which was 
estimated from a comparison of the width obtained using 
NLO and NNLO PDFs 0. Varying M w by the uncer- 
tainty of ±29 MeV from the central value of 80.403 



GeV changes Tw by =f9 MeV in each channel. 

The impact of higher-order EW corrections is deter- 
mined by comparing simulated samples of W — > £vy and 
W — > £vjj events generated by PHOTOS Uncer- 
tainties on T w (r^) of 8 (1) MeV were obtained. The 
correction due to EW box diagrams was determined to 
be 12 MeV in both channels. A systematic uncertainty of 
6 MeV in the box diagram correction was assigned from 
the dependence of the correction on the recoil resolution. 

The uncertainty in simulating lepton identification 
variables was constrained from Z — > £ + £~ decays and re- 
sults in a Tw uncertainty of 10 (6) MeV in the electron 
(muon) channel. Variations in the simulation of the de- 
tector acceptance results in a further small uncertainty 
of 3 (4) MeV in the electron (muon) channel. Table U 
summarizes the sources of uncertainty described above. 

A binned likelihood fit to simulated Mj spectra with 
Tw as a free parameter over the region 90 < Mt < 
200 GeV gives T w = 2118 ± 60 sta t MeV for the elec- 
tron channel and Tw — 1948 ± 67 s t a t MeV for the muon 
channel. Figure 1 shows the Mt distributions of the 
data with the best fits. The electron and muon re- 
sults have a common uncertainty of 27 MeV and are 
combined using the BLUE method [13] to give = 
2032 ± 45 s tat ± 57 sys t MeV. The combination has a x 2 of 
1.6 and a total uncertainty of 73 MeV. No statistically 
significant difference is found between fits using only pos- 
itively and only negatively charged leptons. As a cross- 
check, the W width was also determined from a fit to the 
charged-lepton transverse momentum, which has a dif- 
ferent sensitivity to many of the systematics, and a value 
of Tw consistent with the Mt fit at the < la level was 
obtained. 

The result presented in this Letter is the most precise 
direct measurement of the W width. It can be combined 
with published Tevatron direct width measurements Q 
to give a hadron collider average of = 2056±62 MeV. 
A further combinatioon with the preliminary value ob- 
tained from e+e~ collisions, T w = 2196 ± 84 MeV [23j| . 
gives a new world average value of Tw — 2106 ± 50 MeV, 
in good agreement with the SM prediction. 
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TABLE I: The sources of uncertainty (in MeV) on Tw for the 
W — > ev and W —> [iv measurements. If there is a correlated 
source of error between the two measurements its contribution 
to each measurement is listed in the third column, labeled C. 



Intcrministerial de Ciencia y Tecnologia, Spain; the Eu- 
ropean Community's Human Potential Programme; the 
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Source 




At'"' 
W 


C 


Statistics 


60 


67 




Lepton E or p scale 


21 


17 


12 


Lcpton E or p resolution 


31 


26 




Electron energy loss simulation 


13 






Recoil model 


54 


49 




w 
Pt 


7 


7 


7 


Backgrounds 


32 


33 




PDFs 


20 


20 


20 


M w 


9 


9 


9 


EW radiative corrections 


10 


6 


6 


Lepton ID /acceptance 


10 


7 




Total Syst. 


79 


71 


27 


Total (Stat. + Syst.) 


99 


98 


27 
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r w = 1 948 ± 67 MeV 
X 2 /dof [fit range] =17/21 
X 2 /dof [full range] = 21/29 
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(b) 
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FIG. 1: The transverse mass distributions of the W — > ev 
data (a) and W — > fiv data (b) compared to the best fit. 
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